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Abstract 

A sorting network is a shortest path from 12 • • • n to n ■ ■ ■ 21 in 
the Cayley graph of S n generated by nearest-neighbor swaps. For 
m < n, consider the random m-particle sorting network obtained by 
choosing an n-particle sorting network uniformly at random and then 
observing only the relative order of m particles chosen uniformly at 
random. We prove that the expected number of swaps in location j 
in the subnetwork does not depend on n, and we provide a formula 
for it. Our proof is probabilistic, and involves a Polya urn with non- 
integer numbers of balls. From the case m = 4 we obtain a proof of a 
conjecture of Warrington. Our result is consistent with a conjectural 
limiting law of the subnetwork as n —* oo implied by the great circle 
conjecture Angel, Holroyd, Romik and Virag. 

1 Introduction 

Let S n be the symmetric group of all permutations o = (cr(l), . . . , cx(n)) on 
{1, . . . , n}, with composition given by (ar)(i) := cr(r(z)). For 1 < s < n — 1 
denote the adjacent transposition or swap at location s by r s := ( s s + 1) — 
(1, 2, . . . , s + 1, s, . . . , n) G S n . Denote the identity id := (1, 2, ... , n) and 
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Figure 1: An illustration of the 5-particle sorting network u = 
(2,1,3,4,2,3,4,2,1,2), together with the 3-particle subnetwork uj\a = 
(2, 1, 2) induced by the subset of particles A = {1, 2, 4}. 

the reverse permutation p := (n, . . . , 2, 1). An n-particle sorting network 

is a sequence u = (si, . . . , s^), where N := Q), such that 

T ai T S2 ■ ■ ■ T SN p. 

For 1 < t < N we refer to s t = s t (uj) as the tth swap location, and we call 
the permutation a t = <J t (uj) := r Sl • • -r St the configuration at time t. We 
call o" f _1 (z) the location of particle i at time t. 

Given an n-particle sorting network uj and a subset A of {1, . . . , n} of size 
m, the induced subnetwork uj\a is the m-particle sorting network obtained 
by restricting attention to the particles in A. More precisely, if the elements 
of A are ai < a2 < ■ • ■ < a m , delete from each configuration <j t of uj all 
elements not in A, and replace with i, to give a permutation in S m , then 
remove all duplicates from the resulting sequence of permutations; the result 
is the sequence of configurations of uj\a- See Figure [j]. 

The uniform sorting network uo n is a random sorting network chosen 
according to the uniform measure on the set all n-particle sorting networks. 
For m < n, the random m-out-of-n subnetwork uj 1 ^ is the random m- 
particle sorting network (o; n )|^, where u n is a uniform n-particle sorting 
network, and A is an independent uniformly random m-element subset of 
{1 »}• 

Uniform sorting networks were investigated in || , leading to many strik- 
ing results and conjectures. (A different probability measure on sorting net- 
works was considered in ||.) Our main result is the following surprising 
fact about random subnetworks. We denote the falling factorial (a) r := 
a(a — 1) • • • (a — r + 1) (so r! = (r) r ). 
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Theorem 1. Letm < n. In the random m- out- of -n subnetwork, the expected 
number of swaps in location j does not depend on n, and equals 

E #{t : St «) = j] = °'7 )i ";{^ j ^ ) 1 r H , 1 < j < m - 1. 

(j - l)\{m- j - l)\ 

Given only that the left side does not depend on n, the formula on the 
right side may be recovered by reducing to the case n = m, which gives = 
u> m , and using known results on the uniform sorting network (specifically, 
Proposition ^| below). 

It is natural to seek generalizations of Theorem [I]. For example one might 
ask whether the law of : st(u}^) = j} is the same for each n > m. This 
is true for m = 3 (indeed the law of is the same for all n > 3), but fails 
for m = 4 with n — 4,5. 

From the case m = 4 of Theorem [l| we deduce the following result, which 
was conjectured by Warrington ||. (We abbreviate (si, . . . , sn) to S\ ■ ■ • sn)- 

Corollary 2. For all n > 4, the probability that the random A-out-of-n sub- 
network lies in {123212,321232,212321,232123} is 1/4. 

The present work was triggered by Warrington's conjecture. Corollary 
has a natural interpretation in terms of geometric sorting networks, which we 
define next. Consider a set of n points in M? with no three collinear and no 
two in the same vertical line. The associated geometric sorting network 
is defined as follows. Label the points Xi, . . . ,x n in order of their projections 
onto the horizontal axis. For all but finitely many angles 9, the projections 
onto the line through in direction 9 fall in an order x a h), . . . , x a r n \ corre- 
sponding to a permutation a = ag of the original order, and as 9 is increased 
from to 7r, these permutations form the configurations of a sorting network. 

The four networks listed in Corollary |2| are precisely those geometric 
networks in which one point is in the convex hull of the other three. It 
turns out that all n-particle sorting networks are geometric for n < 4, but 
not for n > 5, as proved in 0. In fact it is proved in that the uniform 
sorting network is non-geometric with probability tending to 1 as n —>■ oo; 
on the other hand, a principal conjecture of is that in a certain sense the 
uniform sorting network is approximately geometric. 

The conjectures in [[| lead to the following precise prediction for the 
limiting law of the random m-out-of-n subnetwork as n —>■ oo. We will prove 
that Theorem |l] is consistent with this conjecture. 



3 



Conjecture 3. Let X±, . . . , X m be independent identically distributed random 
points in M? chosen according to the Archimedes density 

1 

27r \J\ — x 2 — y 2 

on the disc x 2 +y 2 < 1, and let u) m be the associated random geometric sorting 
network. The random m-out-of-n subnetwork satisfies the convergence in 
distribution 

n d ^ 
ui m — ► u m as n — ► oo. 

Conjecture | is implied by 0, Conjecture 3] (and this is implicit in the 
discussion at the end of [§, Section 1] and ||, proof of Theorem 5]). The 
conjecture implies that any statistic of the law of cj™ should converge to the 
appropriate limit; we establish that this indeed holds for the expected value 
in Theorem [I]. 

Proposition 4. Let uj m be the random geometric sorting network of Con- 
jecture |3j. The expected number of swaps in u) m at location j equals the right 
side in Theorem [Z|. 



2 Proof of main result 



We will use the following key properties of uniform sorting networks. 

Proposition 5 ([@]). Consider a uniform n-particle sorting network, and 
write N = (™) . 

(i) The random sequence of swap locations is stationary. That is, 
(si, . . . , sn-i) and (s 2 , • • • , sn) are equal in law. 

(ii) The probability mass function p n of the first swap location is given by 

Proposition [5] is proved in || Theorem l(i) and Proposition 9]. For the 
reader's convenience we also summarize the arguments here. Part (i) follows 
immediately because (si, . . . , sjv) l— ^ ( s 2, • • • , $n, n — Sj) is a permutation of 
the set of all n-particle sorting networks. Part (ii) requires more technology. 
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By a bijection of Edelman and Greene 0, the location of the first swap 
is equal in law to the position along the diagonal of the largest entry in a 
uniformly random standard Young tableau of shape (n — 1, . . . , 2, 1). The 
mass function of the latter may be computed by using the hook formula 
of Frame, Robinson and Thrall || to enumerate Young tableaux with and 
without a given cell on the diagonal. 

The stationarity of the uniform sorting network will play a key role in 
our proof of Theorem |I[ We remark that the random m-out-of-n subnetwork 
u>m is not in general stationary; even ul is a counterexample, as noted in || . 
Stationarity apparently also fails for the random geometric sorting network 
u m of Conjecture [3] (according to simulations), and therefore it presumably 
fails to hold asymptotically for c<j,™ as n — > oo. 

Our proof of Theorem [l] will proceed by relating the mass function p n to 
a Polya urn process, which is defined as follows. An urn contains black and 
white balls in some numbers (which for our purposes need not be integers). 
At each step, one new ball is added to the urn; if the urn currently contains 
w white and b black balls, the next ball to be added is white with probability 
w/(b + w), otherwise black. 

Lemma 6. Consider a Polya urn that initially contains l| black and l| 
white balls. 

(i) The random sequence of colors of added balls is exchangeable (i.e. in- 
variant in law under all permutations affecting finitely many elements). 

(ii) After n — 2 balls have been added, the probability that k — 1 of them are 
white equals p n {k). 

(Property (i) is well known, for arbitrary initial numbers of balls). 

Proof. The probability of adding k — 1 white followed by n — k — 1 black balls 
is 

1§2§ tlly^A ?L n-k- \ _ 2{k- |)fc_i(n - k - |) n _ fc _x 

(1) 

Moreover, the probability of adding k — 1 white and n — k — 1 black balls in 
any given order also equals (|l]), since we obtain the same denominators, and 
the numerators in a different order. This gives the claimed exchangeability. 
Therefore, the probability that k — 1 of the first n — 2 balls are white is (|I]) 
multiplied by ClZi)i which equals the right side in Proposition^ (ii). □ 
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Let /i™ k be the mass function of a hypergeometric distribution, i.e., let 
^mfcW b e ^ ne probability of obtaining i white balls when m are chosen at 
random without replacement from an urn containing n balls of which k are 
white. So 

/,n / -\ _ in /-\ _ [iJ \m-i) 
,l m,k\ l ) — n k,m\ L ) ~ / n \ 

(where as usual we take m = if b [0, a]). 

Lemma 7. For integers n, m,j satisfying m < n and 1 < j < m — 1, 

J2P^ h ™\k-lU - 1) = PmO*)- 

In particular the left side does not depend on n. 

Proof. Consider the Polya urn of Lemma |6|. When n — 2 balls have been 
added, suppose m — 2 balls are chosen at random from these n — 2. Then 
the left side is the probability that j — 1 of those chosen are white. By the 
exchangeability in (i), this probability remains the same if we condition on 
the event that the chosen balls are the first m — 2 to be added to the urn, 
but then the probability is clearly p m (j) by Lemma | (ii). □ 

We remark that a direct computational proof of Lemma [7] is also possible, 
using induction on n. 

Proof of Theorem |||. Consider the random uniform m-out-of-n subnetwork 
u;™ = (u n )\j[. Let q(n,m, k,j,t) be the probability that the ith swap in the 
n-particle network uj n occurs in location k, and that this swap corresponds 
to some swap in location j in the m-out-of-n-network. By the stationarity 
in Proposition |5|(i), and the fact that selecting m items at random from n is 
invariant under permutations of the n items, q is constant in t. On the other 
hand we have 

q(n,m,k,j,l)=p n (k) h n m ^2) h n ~\ k _S - 1), 

since given that the first swap in uj n has location k, the event in question oc- 
curs if and only if the m chosen elements comprising A include the pair 
k,k + 1, and exactly j — 1 of 1, . . . , k — 1. Now the required expecta- 
tion is Yl/k&z St=i l(. n i m i k->3i 0- ^ the above observations, together with 
Lemma [?| and the fact that h 1 ^ 2 (2) = (™) / (™) , this sum equals (^)pm(j)- □ 
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3 Proofs of additional results 



Proof of Corollary [|. It is easy to check that of the 16 4-particle sorting 
networks, the given 4 each have 3 swaps in location 2, while the remaining 
12 each have 2. But Theorem [I] gives that the expected number of swaps in 
location 2 is 9/4 = (1/4)3 + (3/4)2. □ 

Proof of Proposition [|. Let X±, . . . , X m be i.i.d. with Archimedes density as 
in Conjecture |3|. We start by noting two properties. First, the projection of 
Xi onto any fixed direction has uniform distribution on [—1, 1]. Second, the 
signed distance from of the line through X\ and X 2 has semicircle law, i.e. 
density function -yl — r 2 on [—1,1]. (See proof of Theorem 5]). 

Since each pair of particles swaps somewhere in Q m , it suffices to compute 
the probability that a given pair, say those corresponding to X 1; X 2 , swap in 
location j (and then multiply by (™)). This swap occurs when the rotating 
projection line is perpendicular to the line through X\ and X 2 , at which time 
the projections of X\ and X 2 coincide, at a point R with semicircle law. This 
swap is at location j precisely if j — 1 of X 3 , . . . , X m are projected to the left 
(say) of R; but the projections of these points are uniform and independent 
of R. Thus the required expectation is 

Leaving aside multiplicative constants and applying the change of variable 
t = (r + l)/2, the integral reduces to a standard Beta integral (see e.g. |7|, 
p. 148]): 







T(m + 1 



Using r(|) = a/tt, a routine computation then shows that (0) equals the 
right side in Theorem [D, □ 

We remark that the last computation may be viewed as an asymptotic 
version of Lemma [7|, in the limit n — > oo. 
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Open questions 



(i) Does the law of the random m-out-of-n sorting network converge as 
n — > oo, for fixed ml (Conjecture |] makes a specific prediction about 
the limit, but even its existence it is not known.) 

(ii) Our use of the Polya urn can be viewed as a natural way to couple the 
law of the first swap location si(u n ) for different values of n - indeed 
the coupling has the property that si(u n+ i) — si{uj n ) e {0, 1}. Is there 
a natural way to couple the entire uniform sorting networks u n and 
u; n+1 ? For example, can it be done in such a way that u n = (u n+ i)\j3, 
for some random n-element set £>? 
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